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We discuss the role of elastic stress in the statistical properties of elastic turbulence, realized by 
the flow of a polymer solution between two disks. The dynamics of the elastic stress are analogous 
to those of a small scale fast dynamo in magnetohydrodynamics, and to those of the turbulent 
advection of a passive scalar in the Batchelor regime. Both systems are theoretically studied in 
literature, and this analogy is exploited to explain the statistical properties, the flow structure, 
and the scaling observed experimentally. The following features of elastic turbulence are confirmed 
experimentally and presented in this paper: 

(i) the rms of the vorticity (and that of velocity gradients) saturates in the bulk of the elastic 
turbulent flow, leading to the saturation of the elastic stress. 

(ii) the rms of the velocity gradients (and thus the elastic stress) grows linearly with Wi in the 
boundary layer, near the driving disk. The rms of the velocity gradients in the boundary layer is 
one to two orders of magnitude larger than in the bulk. 

(iii) the PDFs of the injected power at either constant angular speed or torque show skewness 
and exponential tails, which both indicate intermittent statistical behavior. Also the PDFs of the 
normalized accelerations, which can be related to the statistics of velocity gradients via the Taylor 
hypothesis, exhibit well-pronounced exponential tails. 

(iv) a new length scale, i.e the thickness of the boundary layer, as measured from the profile of the 
rms of the velocity gradient, is found to be relevant for the boundary layer of the elastic stresses. 
The velocity boundary layer just reflects some of the features of the boundary layer of the elastic 
stresses (rms of the velocity gradients). This measured length scale is much smaller than the vessel 
size. 

(v) the scaling of the structure functions of the vorticity, velocity gradients, and injected power is 
found to be the same as that of a passive scalar advected by an elastic turbulent velocity field. 

PACS numbers: 83.50.-v, 47.27.-i,47.50.+d 
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I. 



INTRODUCTION 



Long polymer molecules added to a fluid make it elastic and capable to store stresses that depend on the history of 
deformation, thereby giving the medium a memory As was recently demonstrated, elastic stress can overcome the 
dissipation caused by relaxating polymer molecules, and can cause complicated and irregular motion in curvilinear 
shear flows 0,0]. The discovery of elastic turbulence, a random flow in dilute polymer solutions at arbitrary low 
Reynolds numbers, Re, calls for quantitative studies of all aspects of the phenomenon, both experimentally and 
theoretically. 

Although the notion of turbulence itself is widely used in scientific and technical literature, there is no commonly 
accepted definition. Turbulent flow is usually identified by its main features. Turbulence implies fluid motion in a 
broad range of temporal and spatial scales, so that many degrees of freedom are excited in the system. There are no 
intrinsic characteristic scales of time and space in the flow, except for those restricting the excited temporal and spatial 
domains from above and below. Turbulent flow is also usually accompanied by a significant increase in momentum 
and mass transfer. That is, flow resistance and rate of mixing in a turbulent flow become much higher than those 
would be in an imaginary laminar flow with the same Reynolds number. 

Random flows of dilute polymer solutions too, cause a sharp growth in flow resistance, exhibit an algebraic decay 
of velocity power spectra over a wide range of scales, and provide a way for effective mixing These properties 

are analogous to those of hydrodynamic turbulence. Such formal similarities give a reason for calling such flows 
"elastic turbulence". However, the similarities do not imply that the physical mechanism that underlies the two kinds 
of random motion is the same. Indeed, in contrast with inertial turbulence at high Re, which occurs due to large 
Reynolds stresses @, large elastic stress is the main source of non-linearity and the cause of elastic turbulence in low 
Re flows of polymer solutions @ . One can suggest that in a random flow driven by elasticity, the elastic stress tensor 
t p should be the object of primary importance and interest, and that it may be appropriate to view elastic turbulence 
as turbulence of the t p field. It would then be more relevant and instructive to explore the spatial structure and 
the temporal distribution of this field, but, unfortunately, currently no experimental technique allows a direct local 
measurement of it in a turbulent flow. On the other hand, properties of the r p field in a boundary layer can be 
inferred from measurements of injected power, whereas its local properties can be evaluated from measurements of 
spatial and temporal distributions of velocity gradients. 

In this regard, we would like to point out the deep similarity between elastic turbulence and the small-scale turbulent 
magneto-dynamo discussed theoretically by @, that will be further clarified in the next Section. Both the small-scale 
magneto-dynamo @] as well as elastic turbulence arise when the flow is three-dimensional, random in time and 
spatially smooth. Analogously to elastic turbulence, the small-scale magneto-dynamo exhibits a viscosity-dominated, 
low Re turbulence that is much faster than the large-scale one. The small-scale turbulent dynamo is produced by the 
random stretching of nearly frozen-in magnetic field lines by the ambient velocity field, which is random and spatially 
smooth, with a growth rate of the order of the inverse viscous eddy-turnover time. The stretching of the magnetic field 
amplifies exponentially the field strength, at the rate of turbulent stretching, and produces dynamically significant 
fields before the large-scale fields can grow appreciably 0|. The advection of the field then spreads the magnetic 
energy over a sub-viscous range of scales and piles it up at the resistive scale. Below that scale, diffusion dominates 
and causes an exponential decay of the energy. This fast dynamo mechanism, due to stretching and folding, generates 
magnetic field lines which remain straight and fairly well-aligned. Similarly, passive scalar mixing in the Batchelor 
regime gives rise to the so-called collinear anomaly Q (see further for details). The fundamental difference between 
the fast small-scale magneto-dynamo and the elastic turbulence lies in what happens when the small-scale magnetic 
field becomes sufficiently strong for the Lorentz force to react back on the flow. The different dissipation mechanism 
in the two systems leads to a different saturation mechanism. In the latter system, saturation is reached when the 
dynamo stretching that produces the feedback on the velocity field and the homogeneous linear damping balance, 
whereas this balance is missing in the former case. 

The microscopic phenomenon responsible for the macroscopic elastic turbulent flow of dilute polymer solutions 
is the so-called coil-stretch transition, discussed originally by Lumley ^3] and then elaborated by de Gennes [Tlj . 
The essence of the coil-stretch transition is as follows: all the polymer molecules are normally, in absence of a flow 
gradient, coiled into compact and nearly spherical tangles with radii much smaller than the total polymer length. 
In the presence of random flow gradients, the polymer molecules are deformed into ellipsoid-like tangles. The mean 
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deformation is moderate as long as the stretching rate of the flow is smaller than the inverse polymer relaxation time. 
When the stretching rate exceeds the inverse polymer relaxation time, most of the molecules are strongly stretched 
to sizes up to the order of the total polymer length, since stretching is limited only by non-linear effects and back 
reaction of the flow. There is an abrupt change in the polymer conformation occurring at large stretching rates, 
called the coil-stretch transition. The transition has remarkable macroscopic consequences. While in the coiled state 
the influence of polymers on the flow is negligible and the polymer solutions can be treated as Newtonian fluids 
(without memory), in the stretched state polymers produce an essential feedback on the flow, and the properties of 
the solution become strongly non-Newtonian. Let us stress that the coil-stretch transition can occur even at small 
Reynolds numbers or at low level of hydrodynamic nonlinearity. The elasticity of the polymers itself (which manifests 
itself above the coil-stretch transition) is a source of nonlinearity. This explains how a chaotic state, which requires a 
high level of nonlinearity, can arise even at low Re in elastic turbulence. 

Another important characteristic of the elastic turbulence regime to be investigated is the scaling of the velocity 
field. The power spectrum of velocity fluctuations measured in the experiments 0,13] shows a power-law behavior f~ s 
with 5 > 3, where / is the frequency. Since the power spectrum decays strongly with the frequency (and thus with 
the wave number) , the fluctuations of the velocity and of the velocity gradients are both determined by the integral 
scale, i.e. the size of the vessel. This explains why elastic turbulence is essentially a smooth random flow, dominated 
by the strong nonlinear interaction of few larg e-scale modes. A theoretical explanation of the power velocity spectrum 
was provided by Fouxon and Lebedev The spectrum k~ s with S > 3 implies that elastic turbulence flow can 

be assimilated to a Batchelor regime flow, i.e. a temporally random and spatially smooth flow 0|. This type of a 
random flow appears in hydrodynamic turbulence below the dissipation scale. 

In a spatially smooth flow, velocity differences are proportional to distances between two fluid particles. As a 
result, particles separate or approach exponentially in time. The statistics of stretching and contraction of the fluid 
element can be described by the formalism of dynamical systems, using the notion of Lyapunov exponents that are 
the exponential rates of divergence. The number of the Lyapunov exponents is equal to space dimensionality, and 
they are ordered in such a way that the first Lyapunov exponent is the largest and describes the separation rate of 
two fluid particles. Based on the Lagrangian approach, the so-called collinear anomaly was predicted (first, within 
the Kraichnan model and then for the general case): a spatially smooth flow locally preserves straight lines This 
anomaly, translated into the spatial distribution of the passive scalar, suggests that parallel strips are the major 
structural elements in the Batchelor regime. The collinear anomaly is also the source of strong intermittency and 
of the absence of scale invariance observed at scales larger than the forcing scale 0. This simple yet fundamental 
feature was observed experimentally for a passive scalar in elastic turbulence 0, Q-Ql- 

Such a Batchelor flow is therefore an efficient mixer due to the exponential divergence of Lagrangian trajectories. 
This property was directl y co nfirmed experimentally in curvilinear channel flows of polymer solutions, where elastic 
turbulence was excited |l,[lj] The passive scalar (dye) mixing had there an exponential character along the channel. 
Additionally, the Lagrangian technique was used (a) to predict the exponential tails of the passive scalar distribution 
function in the steady regime, as was later confirmed experimentally and numerically; (b) to account for a complex 
interplay between Lagrangian stretching and diffusion, thus showing that the PDF of the scalar d issip ation field 
has a stretched-exponential tail [T^ j: (c) to explain strong intermittency in the scalar decay problem |lft Il7j. where 
the decay rate is exponential in time and the exponent saturates to a constant for high-order moments. While the 
experiment 0| confirmed the exponent saturation, it contradicted some other conclusions of the theory developed 
for infinite domains lift Il7l . This observation motivated theorists to develop a new (and much-needed) theory of 
mixing in finite vessels (l£| . This theory explained that long-time properties of scalar decay are determined by the 
stagnation regions of the flow, in particular boundary layers. A nontrivial dependence of the decay rate on diffusivity 
was predicted, as remarkably confirmed by the experiment on mixing in elastic turbulence 0]. These experimental 
discoveries and subsequent theoretical developments show the need to readjust the thinking in fundamental fluid 
mechanics (that mixing of highly viscous fluids is difficult) but also open up a whole new world of practical applications. 

Elastic turbulence gives the possibility to dramatically increase mixing rates (and, correspondingly, rates of chemical 
reactions) of polymer melts in industrial reactors and in small-scale flows (via capillaries, for instance). Indeed, to 
mix fluids effectively one needs irregular motion, which brings into near contact distant portions of the fluids. This 
method of mixing begins to find applications in microfluidics and biophysics ^| . It is likely that Nature itself uses 
the trick of small-scale mixing by elastic turbulence. Biological discoveries and applications along this direction are 
likely to follow. 

In this paper we present experimental results concerning statistics of the injected power and statistics and structure 
of the velocity field in an elastic turbulent flow, in its bulk as well as near a wall. By using the analogy between passive 
scalar advection in the Batchelor regime and advection of elastic stresses in elastic turbulence, we try to elucidate 
the role of the latter, particularly in connection with the intermittency and skewness of the injected power statistics, 
with the statistics and saturation of rms of vorticity (velocity gradients) in the bulk, and with various properties of 
the velocity profile in the elastic turbulent flow regime, such as scaling and properties of the boundary layer. 
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Our experimental program addresses the following questions: 

(i) Does elastic stress saturate in the bulk as predicted theoretically, and what is the saturation value? 

(ii) What is the mechanism for pumping energy from the driven plate into the bulk? 

(iii) What is the probability distribution of the elastic stresses? 

(iv) What is the velocity profile in the cell, and does an analog of a boundary layer exist? 

(v) If the boundary layer exists, which length scale characterizes its thickness, and what is the nature of this scale? 

(vi) Are there similarities in the statistical and scaling behavior with the problem of turbulent advection of a passive 
scalar? 

The paper is organized as follows. We start the next section presenting a hydrodynamic description of dilute polymer 
solution flows and derive non-dimensional parameters to characterize the flows that follow from the equations. The 
variation of one of these control parameters, responsible for the elastic properties of a fluid, can lead to elastic 
instability in the swirling flow between two disks, distinguished by the presence of curvilinear trajectories. The 
stretching of polymer molecules, leading eventually to the coil-stretch transition, is further discussed in section III Bl 
The theoretical criterion for the elastic instability in different flows, together with its experimental verification, are 
discussed in section ITTCl How this instability is responsible for the scenario dubbed elastic turbulence is addressed in 
section lTl Dl The experimental measurement techniques used to characterize the flow are presented in section lTlI Al and 
to complete the basics, the rheometric properties of polymer solutions used are given in section ITlI Bl while previous 
experimental observations of elastic instability in the von Karman flow are summarized in section ITlI CI a nd compared 
to our findings. In section |^ we give a complete account of the results of the measurements. Different aspects of 
the flow are investigated with different techniques: we concentrate on the statistics of global power fluctuations in 
section IIV Al on the topology of the flow below and above transition and its relation with the observed statistics 
in llVBI on profiles of velocity and velocity gradients, and the emergence of boundary layers in section llV CI All the 
new experimental facts are summarized and discussed in in section llV Dl where the role of elastic stress, a recent 
theory of elastic turbulence, and comparative considerations about elastic versus hydrodynamic turbulence are made. 
Finally, in section the results are summarized and conclusions are drawn. 



II. HYDRODYNAMIC AL THEORY OF ELASTIC TURBULENCE 
A. Hydrodynamical description of dilute polymer solution flows 

Solutions of flexible high molecular weight polymers are viscoelastic liquids, and differ from Newtonian fluids in 
many aspects 0]. The most striking elastic property of polymer solutions is, probably, the dependence of mechanical 
stresses on the history of the flow. The stresses do not immediately become zero when the fluid motion stops, but 
rather decay with some characteristic relaxation time, A, which can be well above a second. When a polymer solution 
is sufficiently diluted, its stress tensor r can be divided into two parts, r = t s + t v . The first term, r s , is defined by 
the viscosity of the Newtonian solvent r] s and by the rate of strain in the flow, t s = rj s [V • V + (V • V) T ] . The elastic 
stress tensor, t p , is due to the polymer molecules, which are stretched in the flow, and depends on the history of the 
flow. The equation of motion for a dilute polymer solution has the form Q 



^ + (V.V)V=-^+(^Wv + ^, (1) 
dt P \P J P 

where p is the pressure and p is the density of the fluid. One can see that r p enters the equation of motion linearly. 
The equation has a non-linear term, (V • V)V, which is inertial in its nature. The Reynolds number Re defines 
the ratio of this non-linear term to the viscous dissipative term, ^V 2 V. Therefore the degree of non-linearity of 
the equation of motion can still be defined, even for a polymer solution, by the Reynolds number, Re = VLphu. 
Turbulence in fluids at high Re is a paradigm for strongly nonlinear phenomena in spatially extended systems [fitl2fl|. 

The simplest model incorporating the elastic nature of the polymer stress tensor, r p , is a Maxwell type constitutive 
equation p| with a single relaxation time, A, 



(2) 
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Here Dr p /Dt is the material time derivative of the polymer stress, and rj p = r\ — rj s is the polymer contribution to 
the total viscosity 77. An appropriate expression for the time derivative Dr p /Dt has to take into account that the 
stress is carried by fluid elements, which move, rotate and deform in the flow. The translational motion implies an 
advection term (V- V)t p , while the rotation and deformation of the fluid particles leads to contributions like (V- V)x p 
or t p (V • V) pj. Therefore, some nonlinear terms, in which r p is coupled to V, appear in the constitutive relation 
along with terms linear in them. A simple model equation for Dr p /Dt, which is commonly used for description of 
dilute polymer solutions, is the upper convected time derivative, 



^ = 2Z£ + (V • V)x p - t p (V • V) - (V • Vfr p . (3) 

The equations HIM) together with the expression for t s constitute the Oldroyd-B rheological model of polymer solu- 
tions [j]. One can see that the nonlinear terms in the constitutive equation (Eqs. (03j)) are all of the order X(V/L)t p . 
The ratio of those non-linear terms to the linear relaxation term, r p , that represents dissipation, is given by a dimen- 
sionless expression X(V/L), called the Weissenberg number, Wi, and expresses the ratio of the relaxation time to the 
characteristic flow time. The relaxation term t p is somewhat analogous to the dissipation term in the Navier-Stokes 
equation. 

One can expect mechanical properties of the polymer solutions to become notably nonlinear at sufficiently large 
Weissenbergnumbers. Indeed, quite a few effects originating from the nonlinear polymer stresses have been known for 
a long time pj. In the simple shear flow of a polymer solution, there is a difference between normal stresses along the 
direction of the flow and along the direction of velocity gradient. At low shear rates, the first normal stress difference, 
Ni, is proportional to the square of the shear rate. For curvilinear flow lines, this gives rise to a volume force acting 
on the liquid in the direction of the curvature, called the "hoop stress". Therefore, if a rotating rod is inserted in 
an open vessel with a highly elastic polymer solution, the liquid starts to climb up on the rod, instead of being 
pushed outwards by the centrifugal force. This phenomenon is known as "rod climbing", or "Weissenberg effect" 0|. 
Furthermore, in a purely elongational flow the resistance of a polymer solution depends on the rate of extension in a 
strongly nonlinear way. There is a sharp growth in the elastic stresses when the rate of extension exceeds 1/A, that 
is at Wi > 1. As a result, the ap parent viscosity of a dilute polymer solution in elongational flow can increase by up 
to three orders of magnitude (21| . Both the Weissenberg effect and the growth of the extensional flow resistance were 
most clearly observed in viscous polymer solutions and in flows with quite low Re, when nonlinear inertial effects are 
insignificant 

If there are two potential sources of nonlinearity in the hydrodynamic equations for a dilute polymer solution 
(Eqs. l|H3[l ). inertial and elastic, a natural question arises: how is the hydrodynamical stability of a flow affected by 
each nonlinear contribution? With respect to the values of two non-dimensional parameters, Re and Wi, four cases 
can be considered in viscoelastic hydrodynamics. If both Re and Wi are smaller than unity, the flow is steady and 
laminar. Then in the limit of negligible Wi and sufficiently large Re, the flow loses its stability only because of the 
inertial nonlinearity, which is much larger than the viscous dissipation, and it evolves into the regime of inertial or 
hydrodynamic turbulence of a Newtonian fluid. This remains true as long as the relaxation time A is sufficiently small 
and the value of Wi remains negligible even at the large local velocity gradients appearing in turbulent flow. The 
case of both large Re and Wi is commonly realized in the hydrodynamic turbulence of fluids with small additions of 
long polymer molecules. In such situation the most striking effect observed is that of turbulent drag reduction. The 
last and less discussed case is that of large Wi and negligibly small Re, and is the subject of this paper. 

At first, it is not trivial to realize that some kind of turbulent flow, driven solely by the nonlinear elastic stresses, 
may exist in the absence of any significant inertial effects, at low Re. An important step in this direction was made 
about a decade and a half ago, when purely elastic instabilities were experimentally identified in the rotational flow 
between two plates and in the Couette- Taylor (CT) flow between two cylinders j2^|. 



B. Polymer extension in random flows 

Recently Lumley's theory p~3] was revised, and a quantitative theory of the coil-stretch transition of a polymer 
molecule in a 3D random flow was developed [2J, ■ The dynamics of a polymer molecule are indeed sensitive to the 
motion of the fluid at the dissipation scale, where the velocity field is spatially smooth and random in time. On this 
scale, polymer stretching is determined solely by the velocity gradient tensor, VjV^-, which varies randomly in time 
and space: dtRi = RjVjVi — Ri/X(R) + Q. Here Ri are the end-to-end vector and R is the end-to-end distance for the 
stretched polymer molecule, respectively; X(R) is the polymer relaxation time, which is i?-dependent, and Q is the 
thermal noise. At R <C R ma x, where R ma x is the maximum end-to-end stretched polymer length, the linear regime 
of polymer relaxation is characterized by the polymer relaxation time X re i = A(7 = 0). For molecules significantly 
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more elongated than in the coiled state one can use, e.g., the FENE model with A(i?) = A(7 = 0)(1 — R 2 /Rmax) Ql- 
In a 3D random flow, Y i Y i ^ as a ^ wa y s an eigenvalue with a positive real part, so that there exists a direction with 
a pure elongation flow [26] . The direction and the rate of the elongation flow change randomly, as the fluid element 
rotates and moves along the Lagrangian trajectory. If VjV} remains correlated within finite time intervals, the overall 
statistically averaged stretching of the fluid element will increase exponentially fast in time. The rate of the stretching 
is determined by the maximal Lyapunov exponent, a, of the turbulent flow, which is the average logarithmic rate of 
separation of two initially close trajectories. The value of a is usually of the order of the rms of the fluctuations of 

ol/2 

the velocity gradient, ( J = [g^j 

The stretching of a polymer molecule follows the deformation of the fluid element surrounding it. As a result, the 
statistics of polymer elongations in a random smooth flow depends critically on the value of A • a (or equivalently 

/ i\/ \ rnls i ■ 

A ( g^M ), which plays the role of a local Weissenberg number for a random flow, Wi'. According to the theory (lfl 

l24f . the polymer molecules should become greatly stretched, if the condition A -a > 1 is satisfied. Thus, the coil-stretch 
transition is defined by the relation A cr • a = 1, similarly to that in an elongation flow with the strain rate a [ijj. A 
somewhat surprising conclusion of the theory is that a generic random flow is, in average, an extensional flow in every 
point, with the rate of extension e = a and unlimited Henky strain. A dramatic extension of the flexible polymer 
molecules in the turbulent flow environment, inferred here from bulk measurements of the flow resistance, has been 
recently confirmed by the direct visualization of individual polymer molecules in a random flow ■ 

An appropriate way of describing the coil-stretch transition employs the probability distribution function (PDF) 
of polymer elongations. Below the coil-stretch transition, the distribution is peaked near the equilibrium tangle size, 
determined by thermal fluctuations. Above the coil-stretch transition, the distribution of polymer sizes is changed, 
and becomes peaked at an extension of the order of the total polymer length. The position of the peak is determined 
by an interplay between the stretching rate, the polymer nonlinearity, and the feedback of the polymers on the 
flow. In a random flow, the distribution has an extended (power law) tail related to stretching events that overcome 
polymer relaxation, causing an essential extension of the polymer molecule. The characteristics of the PDF tails and 
their relation to the Lagrangian flow dynamics were theoretically examined in the work by Balkovsky, Fouxon, and 
Lebedev |24[ • According to their recent theory, the tail of PDF of molecular extensions is described by the power law 
P(Ri) ~ R i 13 1 , where (3 ~ (A -1 — a) in the vicinity of the transition. A positive (3 corresponds to the majority of the 
polymer molecules being non-stretched. On the contrary, at (3 < a significant fraction of the molecules is strongly 
stretched, and their size is limited by the feedback reaction of the polymers on the flow [24[ and by the non-linearity 
of the molecular elasticity l2a. Thus, the condition (3 = can be interpreted too as the criterion for the coil-stretch 
transition in turbulent flows [2J|. 

To verify this theoretical picture, experimental studies of polymer stretching in a 3D random flow between two 
rotating disks were conducted in macro- and micro-setups jo,l23- It was experimentally shown that elastic turbulence 
as a macroscopic phenomenon is tied to the coil-stretch transition of polymer molecules that occurs on a microscopic 
level j^jj. In the microscopic experiment j^, the size distribution of DNA molecules in the elastic turbulent flow 
was directly measured by optical visualization methods. The validity of the criterion for the coil-stretch transition 
was quantitatively verified in the experiment by the analysis of both the velocity field in the Lagrangian frame (via 
statistics of the Lyapunov exponents, which defines the average a), and of statistics of single molecule stretching in 
the same flow [27| . A good agreement between theory and experiment was found, both concerning the shape of the 
polymer size distribution and the functional form of its tails. 



C. Elastic instability and its experimental investigation 

During the past decade, the purely elastic instabilities in viscoelastic fluids have been the subject of many theoretical 
and experimental studies, which are partially reviewed in Ref. [28L l2fl| . Purely elastic instabilities of viscoelastic 
flows occur at Wi of order unity and vanishingly small Re. As a result of the instability, secondary vortex flows are 
developed p3j |. and the flow resistance increases |23|. The analysis confirmed that the nonlinear mechanical properties 
of the polymer solution can indeed lead to a flow instability, for which a simple mechanism was proposed. After the 
pioneering work by Larson, Muller and Shaqfeh |23t l-3fl|. purely elastic instabilities were also found in other shear 
flows with curvilinear streamlines. Those included the flow between a rotating cone and a plate and the Taylor-Dean 
flow [i3,[2<|. The original mechanism proposed in Ref. [13] was verified experimentally in Ref. j2j|. The original 
theoretical analysis of Ref. was refined, and more elaborate experiments were carried out. A few new mechanisms 
of flow instability driven by nonlinear elastic stresses were suggested for those geometries. 

The most thorough and detailed studies were conducted on the Couette- Taylor (CT) flow between two coaxial 
cylinders. In spite of the fact that instabilities in viscoelastic fluids were studied for decades, the purely elastic 
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instability in CT flow was first investigated both experimentally and theoretically rather recently |23L l3(i| . The 
mechanism of the elastic instability in the CT system suggested in Ref. [13] is based on the Oldroyd-B model. The 
primary flow in the CT system with the inner cylinder rotating (Couette flow) is a pure shear flow in the r#-plane 
that generates a normal stress difference Ni = Tgg — t„ = 2rj p \( A f r g) 2 and a radial force N\/r per unit volume. Here 
r, 9 are the cylindrical coordinates in the plane perpendicular to the cylinder axis, Tgg and r rr are the components of 
the polymer stress tensor t p , and ^ r g is the only non-zero component of the rate-of-strain tensor, 7. A secondary flow 
in the CT flow includes regions of elongational flow in the r-direction with i = dv r /dr ^ 0. The radial extensional 
flow stretches the polymer macromolecular coils in the r-direction, though it is a small perturbation on top of the 
azimuthal stretching produced by the primary shear flow. However, when stretched in the radial direction, the 
macromolecular coils become more susceptible to the basic shear flow. The coupling between the radial and shear 
flow leads to further increase of shear stresses. Thus, it results in further stretching of the polymer in the ^-direction 
that generates an additional normal stress difference. So, the elastic instability mechanism is based on the coupling 
between the perturbative radial elongation flow and the strong azimuthal shear flow, which results in a radial force. 
The latter reinforces the radial flow. As pointed out in Ref. |2jJ, this transition can be only a finite amplitude (first 
order) transition. The corresponding criterion of the instability is 

K = ^ 4-Wi 2 = 0(1). (4) 

Here d and R\ are the gap and the inner cylinder radius, respectively. The Weissenberg number is defined here 
as Wi = XVlR/d, where ft is angu lar velocity of the rotating inner cylinder. (It was termed as the Deborah number 
in some of the original texts [3CL l3ll The elastic instability occurs when the parameter K exceeds a certain 

threshold value This criterion is valid at sufficiently small values of the polymer viscosity ratio, r) p /ri s , and in 
the limit of small gap ratio d/R\. A more general expression is given in Ref. When both ratios are fixed, the 

elastic instability is defined only by the critical Weissenberg number Wi c - 

As was suggested in Ref. J22(, where the CT flow was discussed, there is some analogy between flow transitions 
driven by elasticity and by inertia. The inertially driven Taylor instability occurs at constant Taylor number 0, [2H, 
Ta = -^-i?e 2 , while the elastic instability is controlled by the parameter K of Eq. J3J For the viscoelastic case, 

the Weissenberg number appears to be analogous to the Reynolds number. The geometric parameter determining 
curvature, d/R\, enters the expressions for both Ta and K. Scales of time and velocity for the purely elastic flow 
transition are given by A and d/X, respectively. As it was shown in Ref. they are analogous to t v d, the viscous 
diffusion time defined as t v d = d 2 jv, and to d/t v d, which provide characteristic scales of time and velocity for the 
inertially driven flow transitions. 

Beyond the analogy, there are important differences between the flow transitions driven by inertia and by elasticity. 
It is an experimental fact that any laminar flow of a Newtonian fluid becomes unstable at sufficiently high Re, 
and all high Reynolds number flows are turbulent. That includes rectilinear shear flows, such as the Poiseuille flow 
in a circular pipe, and the plane Couette flow, which are known to be linearly stable at any Re. In contrast to 
it, purely elastic flow instabilities in shear flows have been observed so far only in systems with curvilinear stream 
lines. All these instabilities are supposed to be driven by the hoop stress, which originates from the normal stress 
differences. Recently, the authors of Ref. presented theoretical arguments and numerical evidence for nonlinear 
instability in viscoelastic plane Couette flow. This theory is waiting its experimental verification and still remains 
quite controversial. 

The difference between the inertial and elastic instabilities originates, of course, from the distinct governing equa- 
tions. There are, however, some purely practical factors that can explain rather well the lack of observations of 
purely elastic flow transitions in rectilinear shear flows. Inertial instabilities in rectilinear shear flows of Newtonian 
fluids occur at quite high Reynolds numbers. These are typically about two orders of magnitude higher than Re at 
which curvilinear shear flows with large gap ratios become unstable. A priori, one may suggest that rectilinear and 
curvilinear shear flows would have a similar relation between Wi at thresholds of the purely elastic flow instabilities 
as well. The problem is that, while it is rather easy to generate high Re flows with Newtonian fluids of low viscosity, it 
is usually impossible to reach the corresponding high values of Wi in shear flows of elastic polymer solutions. That is, 
there are always rather severe practical limitations restricting nonlinearity in elastic polymer stresses in shear flows, 
such as rapid polymer degradation at high Wi. 

A purely elastic instability, apparently due to a similar mechanism, also occurs in other rotational shear flows, 
namely concentric cone-and-plate and plate-and-plate geometries, which are common geometrical configurations for 
rheological tests. In these cases the secondary flow is also driven by a "hoop stress" and by the coupling between the 
primary shear flow and secondary elongation flow in the axial direction. Historically, a transition in the rheological 
measurement was first observed experimentally in these geometries [22], and Magda and Larson suggested that it 
had to be ascribed to an elastic flow instability. Complete sets of numerical and experimental studies of the elastic 



8 



instability in these flows for small aspect ratio vessels were carried on in Ref. [aHHHl- Due to the more complicated 
flow structure than in CT flow, no analytic expression for the stability criterion is available in either case. However, 
some features common for all rotational shear flows can be pointed out. For all systems with a small aspect ratio, the 
threshold of the elastic instability depends on the aspect ratio that is defined differently in each case. The instability 
criterion also depends on the viscosity ratio, though the functional form of this dependence was found only in the CT 
flow. In all systems, the most unstable mode is non-axisymmetric and oscillatory and the transition is discontinuous 
or of the first order (inverse bifurcation) . 

D. Theory of the elastic turbulence and role of the elastic stress 

Let us consider now the limiting case of Wi 3> 1 and Re <C 1 . Along with the apparent similarity in phenomenology 
of both elastic and inertial turbulence, there are also many very important distinctions. One similarity is reflected 
by the fact that the Weissenberg number, which acts as the control parameter in elastic turbulence, depends on the 
system size and on the fluid viscosity, like the Reynolds number in hydrodynamic turbulence. Indeed, the scaling of 
the main characteristics of the elastic turbulence with the system size, fluid velocity, and viscosity is distinctly different 
from hydrodynamic turbulence. For instance, it has been shown that with enough elasticity, one can excite turbulent 
motion at arbitrary low velocities and in arbitrarily small vessels 0,0]. An obvious reason for the difference in scalings 
in the two scenarios is that the physical mechanisms which underlie the two kinds of turbulent motion are themselves 
different. As it is well known, the high flow resistance in inertial turbulence is due to large Reynolds stresses. The 
Reynolds stress tensor is defined as the average value of p(ViVj), where Vi and Vj are different components of the 
flow velocity. In the case of elastic turbulence, the Reynolds stresses are quite small, since Re is low. It follows that 
the high flow resistance in elastic turbulence can only be due to the large elastic stresses, r p |fj. 

As we have already pointed out, the statistical properties of a turbulent elastic flow and the increase in the 
contribution of the elastic stresses to the flow resistance, are associated with the significant polymer stretching in 
the random flow. To support this theoretical idea, some experimental studies of polymer stretching in a 3D random 
flow between two plates were conducted in macro- and micro-setups 0, ■ The first experiment was performed in 
a swirling flow set-up with the high aspect ratio R/d = 1, where R = 30 mm, using a very viscous solvent with 
77 s = 1.36 Pas to suppress inertia jfj. By experimental analysis and estimates of the contributions of the Reynolds, 
viscous, and elastic stresses to the shear stress at the upper plate, we found that the Reynolds stress accounts for 
less than 0.1% of the whole stresses, and that the viscous stress is rather constant. So, as a result of a secondary 
random 3D flow superimposed on an applied primary shear flow between two plates, the polymer contribution to the 
shear stress increases as much as 170 times. If one assumes a linear elasticity of the flow-stretched polymer molecules 
(polyacrylamide, in this case), then the elastic stress causes an extension by 13 times @. 

The molecular theory of polymer dynamics tells that t„ is proportional to the polymer concentration, n, and to the 
average polymer conformation tensor, r Pj y ~ n(RiRj) [y. Then, a significant increase in the contribution of t p to 
the flow resistance in the elastic turbulence is associated, on a microscopic level, with significant polymer stretching 
in a random flow flfl l2i l2fi| . 

The crucial theoretical step towards the description of elastic turbulence relates the dynamics of r p to the dynamics 
of a vector field with a linear damping f\ fl^L l37l | . The elastic stress tensor can be written as uniaxial, i.e. Tik = BiBk, 
if the contribution of thermal fluctuations and polymer nonlinearity to t p can be neglected 0|. Then, one can write 
an equation for Bi in a form which is similar to the equation for magnetic field in magneto-hydrodynamics (MHD) 0|. 
Thus, in the case of a viscoelastic flow, one gets: 

d t B + (V • V)B = (B-V)V - B/A. (5) 

The difference with MHD lies only in the relaxation term that replaces the diffusion term. As already mentioned in 
the Introduction, elastic turbulence is analogous to a fast small-scale viscosity-dominated magneto-dynamo, which 
realizes a sort of low Reynolds number Batchelor regime of turbulence. The latter results from the random stretching 
of the (nearly) frozen-in magnetic field lines by an advecting random flow. Numerical simulations of the former show 
features very similar to those observed in elastic turbulence @. Eq. (0, complemented by the equation of motion, 
written for Re <C 1 as 

VP = p(B-V)B + ?/V 2 V, (6) 

and by appropriate boundary conditions, reveals the elastic instability at Wi = Wi c > 1 USUI, where Wi = 7 A 
and 7 is the shear rate. At Wi > Wi c the instability eventually results in chaotic, statistically stationary dynamics. 
According to Ref. 0], a statistically steady state occurs due to the back reaction of stretched polymers (or the elastic 
stresses in Eq.©) on the velocity field, which leads to a saturation of r p even for linear relaxation. On the other 
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hand, the velocity gradients become smaller with decreasing scale. This means that large-scale fluctuations (of the 
vessel size) are dominant in the flow, and the dynamics are determined by a non-linear interaction of modes on scales 
of the order of the system size. Assuming that viscous and relaxation dissipative terms in Eqs. (JHJl and © are of 
the same order, one gets that the velocity gradients on the small scale, being fixed by the stationarity conditions for 
statistics, are of the order of A -1 So, the elastic stress can be estimated as 

r p = B 2 ~ (i,/p)ViVj ~ (7) 

The fluctuating velocity vector and the stress tensor field (or B) can both be decomposed into large, V and B , and 
small, v and B', scale components. The analysis of the equations for the small-scale fluctuations of both fields leads to a 
power-like decaying spectrum for the elastic stress, which in spherical representation looks like {B^B'j) ~ F(k) ~ fc~ 5 , 
where 5 > 3. It is clear that the field B[ (and so field B^B'j — ViVj) is the passive field in the problem. The 
mechanism of stretching and folding of the elastic stress field by a random advecting flow, with the key dynamo effect 
and homogeneous attenuation of the stress field [23, is rather general. It leads to the power-law spectrum for small- 
scale fluctuations of B' in a chaotic flow, and is directly related to the classical Batchelor regime of mixing j2^|, where 
a passive scalar is advected by a flow which is spatially smooth and random in time. A linear relation between the 
small-scale fluctuations of the fields B' and v allows one to establish the power-law spectrum of for the velocity too, 
which, in spherical representation looks again like E(k) ~ k~ s , where 5 > 3. This agrees well with the experimental 
values S = 3.3 — 3.6 @> H HH • Since the velocity spectrum decays faster than k~ 3 , the powers of both the velocity 
and of the velocity gradients are determined essentially by the vessel size; therefore elastic turbulent flow is spatially 
smooth, strongly correlated on the global scale, and random in time. That is also the main feature of the Batchelor 
regime, in which smoothness in space and randomness in time flow are observed below the dissipative scale [26] . 

We would like to emphasize that hydrodynamic turbulence instead exhibits spatially rough flow on all scales above 
the viscous cut-off, and cannot be approximated by a linear velocity field. Moreover, in contrast to hydrodynamic 
inertial turbulence, the algebraic decay of the power spectrum in elastic turbulence is not related to the energy cascade 
nor to any conservation law, since the main energy dissipation occurs at the largest scales. The mechanism leading 
to the algebraic power spectra for the elastic stress and the velocity in this case is related to the advection and linear 
decay accompanied by stretching of the fluid element carrying the elastic stress. The fast decay of the fluctuation 
power spectrum with k implies a velocity field, where the main contribution to deformation and stirring (stretching 
and folding) at all scales comes from a randomly fluctuating flow field at the largest scale of the system. The suggested 
mechanism for the generation of small scale (high k) fluctuations in the elastic stress is the advection of the fluid 
(which carries the stress) by this fluctuating large scale velocity field. The theory considers the elastic stress tensor 
being passively advected in a random velocity field, just as the passively advected vector B' in the magnetic dynamo 
theory Q. 

The smoothness of the velocity field in elastic turbulence was experimentally tested by investigating the shape of 
the cross-correlation functions of the velocity field 0] . It was found there that the second order spatial derivative 
of the velocity field is roughly an order of magnitude smaller than the first one. Thus, a linear velocity field (locally 
uniform rate of strain) is a good local approximation for the velocity field in elastic turbulence [lflj . 

III. EXPERIMENTAL 

A. Experimental set-up and procedure. 

The experiments were conducted in the apparatus schematically shown in FigQ] It consists of a stationary cylindrical 
container C, with radius R c , mounted on the basis of a commercial rheometer (AR-1000 from TA Instruments) and 
a rotating disk D, with radius Rd, mounted concentrically on the shaft of the rheometer. In order to check the 
sensitivity of the results with respect to the geometrical aspect ratio and to explore a broad range of Weissenberg 
numbers, two versions of this setup were used: the first had R c = 2.2cm, Rd = 2cm and height d — lcm, and the 
second had R c = 4.9cm, Rd — 4.7cm and d = lcm. The two setups will be further referred to as setup 1 and setup 
2. The cylinder C was filled with fluid just up to the disk, whose lower side only was wet. In order to provide 
thermal stabilization, the whole rheometer was placed in a thermally isolated box with through flow of temperature 
controlled air. The rheometer mechanism allowed precise control (within 0.5%) of the angular velocity ft of the disk, 
and simultaneous measurements of the applied torque T (fi-forcing mode), or conversely - control (within 0.5%) of T 
and measurements of Q (T-forcing mode). The average shear stress at the upper plate, t w , was calculated using the 
equation T = t w J rdS, that gave t w = 3T/(2irR d ). (The integration is carried over the upper plate surface). 
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Figure 1: Overview of the swirling flow apparatus: C-fluid container, D-rotating disk, AR-1000-rheometer, CO-cylindrical 
optics, L-laser, M-plane mirror, CCD-digital camera, P-LDV probe, BSA-burst spectrum analyzer, PMT-photomultiplier. 

The moment of inertia of the shaft of the rheometer was I s « lAfiN ■ m ■ s 2 and that of the upper plate, Id, was 
about 61 fj,N • m • s 2 for setup 1, and about 84^iV • m • s 2 for setup 2. The accuracy of the angular speed measurements 
in constant torque mode is about 2% and the accuracy of the torque measurements in the constant speed mode is 
about 1%. One has to point out here that the small rate of the fluctuations of the angular velocity is not a sufficient 
criterion to have a constant speed forcing. In the 17-mode, for instance, (I s + Id)^ should also be much smaller than 
the typical values of torque, T. 

To obtain a detailed characterization of both the spatial structure and the temporal evolution of the flow-fields, three 
different experimental techniques have been alternatively used: particle image velocimetry (PIV), particle tracking 
velocimetry (PTV) and laser Doppler velocimetry (LDV) . The cup C was machined from a single piece of optically 
clear perspex, to allow side and bottom view. The cup had a circular inner and a square outer horizontal cross-section 
(i.e. flat sides) in order to insure a distortion free illumination, and was attached to the rheometer base in such a way 
that view from below was possible through a mirror M, tilted by 45°. As flow tracers, we have used 10/xm fluorescent 
particles for the PIV and PTV measurements, and l/im latex beads for the LDV measurements. For PIV and PTV 
the system was illuminated laterally by a thin (30^m in the center of the set-up and about 130//m at the edges of 
the set-up) laser sheet through the transparent walls of the fluid container at the middle distance between the plates. 
The laser sheet was generated by passing a laser beam delivered by a 300mW argon-ion laser, L, through a block of 
two crossed cylindrical lenses, CO, mounted in a telescopic arrangement. The mirror M was used to relay images of 
the flow to a CCD camera (PixelFly digital camera from PCO). The camera was equipped with a regular video lens 
and mounted horizontally near the rheometer (FigQJ. Images were acquired with 12 bit quantization and resolution 
of 640 x 512 pixels at up to 25 frame/sec and 1280 x 1024 at up to 12.5 frame/sec. 

The main tool for the investigation of the flow was the digital PIV technique. We acquired series of up to 2000 
pairs of flow images with the camera. A custom development of the camera control software allowed us to adjust the 
time delay between consecutive PIV images in relation to the local flow velocity (in order to keep the mean particle 
displacement in the range 5 — 15 pixels). Corresponding to low values of the angular velocity of the upper disk (up to 
about OArad/s) the time delay was 540ms and then, gradually decreased down to 40ms. The total data acquisition 
time was always longer than the Eulerian correlation times of the velocity (about 24 times in the fully developed 
random regime). Time series of velocity fields were obtained by a multi-pass PIV algorithm (T3 [. The accuracy of 
the method was carefully checked by running test experiments with the solvent, in the same range of mean particle 
displacements and at similar illumination conditions. Although the instrumental error increases more or less linearly 
with the mean particle displacement, it never exceeded 5% of the mean displacement. The spatial resolution of the 
PIV velocity measurements was 0.57mm in setup 1 and 0.89mm in setup 2. By post-processing the velocity fields, 
we obtained the profiles of the velocity components, fields of fluctuations of each velocity component, spatial spectra 
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of the velocity fluctuations, velocity gradients and their fluctuations, structure functions of gradients, and Eulerian 
velocity correlation functions. The space-time measurements together with simultaneous global measurements of the 
flow resistance provided a rather complete description of the different flow regimes as a function of Wi (IE ED] ■ 

We also examined Lagrangian properties of the flow In order to check the consistency of the PIV approach 
(particularly when Lagrangian trajectories were measured), flow fields have been alternatively investigated with the 
PTV technique. The first step of the PTV approach is the accurate particle identification in each flow image, by 
correlation with prototypical particle shapes, which can be either extracted from flow images, or defined manually. 
Then, trajectories are reconstructed by joining successive positions of the same particle on subsequent images. This 
has been done by a search algorithm based on an initial "guess" of the mean flow line. 

Although the PIV and PTV tools are very suitable for the investigation of the spatial properties of the flow, 
they provide a limited time resolution and statistics. In order to overcome this limitation, the LDV technique has 
alternatively been used. With the LDV technique, time series of the fluid velocity are measured in a relatively small 
volume (typically 100/i?7i x 50/ito x 50fim). The setup allowed measurements of one component of the velocity by 
LDV with two crossing and frequency shifted beams. By appropriate positioning and orientation of the beam crossing 
region, an azimuthal (longitudinal) velocity component, Ve, could be measured at different r and z. Here (r,6,z) are 
the standard cylindrical coordinates. A commercial LDV system from Dantec Dynamics Inc. was used. Two laser 
beams delivered by a probe (P in Fig.QJ are aligned to cross each other at an angle of 9 = 10° at some point in the 
flow. The light scattered in the backward direction by the flow tracers present in the flow is collected by the probe P 
and delivered to a photomultiplier, PMT, via a fiber optic light guide. The PMT signal is processed by a real time 
burst spectrum analyzer, BSA. The local flow velocity is proportional to the Doppler shift of the signal. The error 
validated velocity data are stored on a computer via a GPIB line. 

Our experiments need to be classified according to meaningful values of the non-dimensional control parameters Re 
and Wi, which depend among the rest on a global shear rate. In a swirling flow between two plates, the shear rate 
is quite inhomogeneous over the fluid bulk, even when the flow is laminar. Therefore, the choice of a representative 
shear rate becomes somewhat arbitrary. We decided to consider the simple expression Q,R/d as a characteristic shear 
rate, and to define the Weissenberg number as Wi = XQR/d. The Reynolds number was defined as Re = QRdp/rj. 

B. Rheometric properties of polymer solutions 

The polymer used was polyacrylamide, PAAm (from Polysciences Inc.), with a molecular weight M w = 1.8 x 10 7 Da, 
as used in Ref. Q,E|Q- Our preparation recipe was as follows: first we dissolved 0.9 g of PAAm powder and 3 g of 
NaCl into 275 ml of deionized water by gentle shaking. NaCl was added to fix the ionic contents. Next the solution was 
mixed for 3 hours in a commercial mixer with a propeller at a moderate speed. The rationale for this step is to cause 
a controlled mechanical degradation of the longest PAAm molecules, in order to "cut" the tail of the molecular weight 
distribution of the broadly dispersed PAAm sample. In a solution with a broad distribution of polymer molecular 
weights the heaviest molecules, which are most vulnerable to mechanical degradation, bring the major contribution 
to the solution elasticity, but may break in the course of the experiment. This can lead to inconsistency of the 
experimental results. We found empirically that the procedure of pre-degradation in the mixer leads to substantial 
reduction of degradation during the experiments and to substantial improvement of their consistency Q . Finally, 9 
g of isopropanol was added to the solution (to preserve it from aging) and water was added up to 300 g. The final 
concentrations of PAAm, NaCl and isopropanol in the stock solution were 3000 ppm, 1% and 3%, respectively. This 
master solution was used to prepare 80ppm PAAm solution in a Newtonian solvent. The Newtonian solvent was about 
65% saccharose in water 

The rheological properties of the solvent and the polymer solution were measured with two different rheometers: 
the AR-1000 from TA Instruments around which the complete experimental setup was built, and a Vilastic 3 from 
Vilastic Scientific. The viscosity of the solvent was found to be r] s = HAmPa ■ s at 22°C, and that of the solution 
rj = 138mPa ■ s at a shear rate of 2s" 1 . The polymer relaxation time A was measured in oscillatory tests at different 
shear rates, with 7 ranging from OAs' 1 to 3.6s -1 . In the limit of 7 — > the relaxation time was A(0) = 4.7 s. The 
solution showed a clear "shear thinning" as a function of the shear rate with the scaling A ~ j~ s , where S ~ 0.3, 
similar to what was found earlier for polymer solutions with lower molecular weight PAAm samples (3ll |. 

C. Experimental observation of elastic instability in von Karman swirling flow 

The functional dependence of the elastic instability threshold on the aspect ratio, e = d/R c , in a flow between two 
plates was first obtained in the experiment by Magda and Larson [22|. They found an inverse dependence of the 
critical shear rate %, or of the critical Weissenberg number, Wi c = 7 C A, on the aspect ratio, i.e. Wi c cx e _1 . The 
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experiment was conducted in the range of low values of e between 0.024 and 0.16. In the later experiments similar 
results were obtained in the same range of e by McKinley and collaborators j3 HE H3| ■ These results were found in 
reasonable agreement with numerical calculations. However, these numerical results presented in Wi c — e coordinates 
can be fitted by Wi c oc e -0- 470 ±0- 006 m a wide range of e between 10 -3 and 1, in odds with the the above mentioned 
dependence Wi c oc e _1 . 

We have investigated experimentally the dependence of the onset of the elastic instability on the aspect ratio of 
our setup. By modifying both the radius of the container, i? c , and the distance between plates, d, the aspect ratio e 
was varied by about 10 times for a rather wide range of e, between 0.091 and 0.882. The onset values of the angular 
velocity, O c , were obtained from the dependence of the averaged torque on the angular velocity of the disk (Fig. EJ. 
By linear approximation of this dependence on both sides of the transition point, the latter can be determined within 
reasonable error bars (see Fig- EK) - The Reynolds number corresponding to the onset of the elastic instability, Re c , 
was of the order of unity in both setups, suggesting that the transition is of purely elastic nature. The results of these 
measurements together with other parameters are presented in Table [I] 
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Table I: Onset values of the relevant physical quantities corresponding to the primary elastic instability 

We remark that the values of the critical Weissenberg number as a function of the as pec t ratio obtained from our 
experimental results agree fairly well with the linear stability analysis presented in [34L l35l l3rl| . as shown in FigEJa. 
On the other hand, the fit to the theoretical functional dependence in a wide range of e between 10 -3 and 1, namely 
Wi c oc e -0 - 47 , is quite different from that obtained from the experimental data: Wi c oc e -°- 79 ± 016 . This discrepancy 
can be attributed to the larger aspect ratios experimentally studied by us. 
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Figure 2: Dependence of the onset of the primary elastic instability on the geometric aspect ratio: (a) Torque as a function 
of angular velocity for the aspect ratio e = 0.88. The full lines are linear fits. The arrow indicates the onset of the primary 
elastic instability, (b) Critical Weissenberg number for the onset of the elastic instability as a function of the aspect ratio. The 
symbols are (color online): (black) squares-i? c = 4.7cm, (red) circles-i? c = 2.2cm, (yellow) triangles-i? c = 1.7cm. The full line 
is a linear fit to the data, Wi c = 0.51e~ - 79 . The dashed line is the fit to the numerical result from 0], Wi c = 0.62 + 5.42e~ ' 47 . 
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Figure 3: (a) Dependence of the scaled average power, P/Pi am , on the control parameter Wi: squares-increasing Wi, circles- 
decreasing Wi. The full line is a guide for the eye, with slope Wi ' 49 . (b) Dependence of the scaled rms of the power fluctuations, 
prms J prms Qn ^. S q Uares _ mcreasm g Wi, circles-decreasing Wi. The full line is a power law, Wi 3 ' 2 . Data were collected in 
set-up 2 in forcing mode. 



IV. PROPERTIES OF ELASTIC TURBULENCE 
A. Power fluctuations and statistics of the injected power 

One of the main features of the transition to elastic turbulence is the substantial growth of the flow resistance 
above the onset of the instability 0,0]. In the case of a von Karman swirling flow, a measure of the flow resistance 
is either the power Pq needed to spin the upper disk at constant angular speed Cl, or the power Pt needed to spin 
the upper disk with constant torque T applied to the shaft of the rheometer. The injected power is defined in any 
case as P = T ■ fi. In our experiments, due to the smallness of Re < 16, the inertial contribution was always low. 
The dependencies of the reduced average injected power, P/Pi am and of the rms of its fluctuations, P rms / ! Pf™ ', in a 
wide range of the control parameter Wi in the ft mode are presented in Fig. El Here P\ am and P/ a ™ s are the injected 
power and its fluctuations before the onset of the elastic instability. 

The data presented in Fig. reveal three distinct flow regimes. For low values of the control parameter Wi, the 
reduced injected power is equal to unity. At Wi ~ Wi c , a primary elastic instability occurs, resulting in a sharp 
increase of the injected power and of the rms of its fluctuations. A further increase of the control parameter causes 
the flow to evolve towards a scaling regime, where both the reduced average and rms of fluctuations of the injected 
power scale algebraically with Wi (Fig.©: P/Pi am oc v^i a49±0 05 and P rms /P[™ oc Wi 3 - 2±0 3 . 

The fluctuations of the injected power were measured for different Wi in the elastic turbulence regime in two modes: 
fi-forcing is presented in Fig. and T-forcing is presented in Fig.0 For each value of Wi the statistics of the power 
fluctuations was collected on 180000 data points evenly sampled in time (At « 38ms). The time series of the injected 
power at different Wi for both modes are presented in both figures. In order to avoid mechanical degradation of the 
the polymer solution during the long data acquisition times, separate experiments were conducted for each value of 
the control parameter. The local (in time) average values of the injected power did not change significantly during the 
total data acquisition times, suggesting that no major degradation occurred. For low values of the control parameter 
in the laminar regime, (Wi < Wi c ), the power fluctuations are only due to the instrumental noise. In the elastic 
turbulence regime, the PDFs of the power fluctuations strongly deviate from Gaussian distributions for both forcing 
modes. The PDFs in the O-forcing mode have a left side skewness, while in the T-forcing mode -a right side skewness. 
Since the fluctuations of the injected power reflect fluctuations of the elastic stress averaged over the upper plate, one 
can conclude that the statistics of the elastic stress too must strongly deviate from Gaussian, which is a signature 
of intermittency. The PDFs of the reduced power fluctuations collapse on a single curve when normalized by their 
maximum value (see Fig.©. Both the dependency on Wi of P / P rms (Fig.©, and the skewness of the PDF in elastic 
turbulence suggest an analogy with similar behavior in hydrodynamic turbulence, though the reasons for the effects 
are entirely different [421.I43I]. 

We further examined the statistics of rare and short time scale events (spikes) in the injected power series; namely 
we tried to quantify the characteristic width of the spikes and the characteristic time between them. To carry out this 
analysis, the following procedure was applied. First, all local minima and maxima were detected in each time series. 
Then only those minima (in the case of constant ft forcing) and those maxima (in the case of constant T forcing) that 
contribute to the skewness of the corresponding PDFs, are selected, i.e. those peaks which are larger than 2P rms . 
The values of the average time t s between two spikes, normalized by the relaxation time, are presented in Fig. 0a), 
as function of Wi in both ft and T modes. 
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Figure 4: (a) Time series (partially shown) of the injected power in constant SI forcing regime, for different Wi. The labels 
are: (l)-Wi = 5, (2)-Wi = 19, (3)-Wi = 24, (4)-Wi = 31.5, (5)-Wi = 40. (b) PDF's of the power fluctuations inn-forcing 
regime for different Wi. The symbols are (color online): squares (orange)- Wi = 5, left triangles (black)- VKi = 19, up triangles 
(red) -Wi = 24, diamonds (green)-Wi = 31.5, circles (b\ue)-Wi = 40. Data were collected in setup 2. 
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Figure 5: (a) Time series (partially shown) of the injected power in constant T forcing regime for different Wi. The labels 
are: (l)-Wi = 5, (2)-Wi = 17, (3)-Wi = 25, (4)-Wi = 34, (5)-Wi = 46. (b) PDF's of the power fluctuations in T-forcing 
regime for different Wi. The symbols (color online) are: down triangles (blue^Wi = 5, up triangles (green)-W / i = 17, circles 
(xed)-Wi = 25, diamonds (orange)-VKi = 34, squares (black^VKi = 46. Data were collected in setup 2. 



We performed a similar analysis on the time series of the local azimuthal velocity measured by the LDV technique, 
such as those shown in the following in Fig. El The values of the the average time between two spikes, normalized 
by the relaxation time, are presented in Fig.EJb) as function of Wi. We observe that the spikes in the local velocity 
signal occur about three time less frequently than those in the injected power series. This results in poor statistics 
and larger scatter, particularly for low values of Wi, but on the other hand is consistent with the fact that the rare 
events in the power time series are actually a result of space averaging (over the entire volume of the flow cell) of 
individual rare events in local velocity time series. 

We also computed the mean width of the spikes determined by the procedure described above. The central part 
of each spike was fitted by a Gaussian function, and all the fit widths c for a given Wi were averaged together. The 
results for the spikes in the series of injected power in the fi-mode are shown in Fig. |H| 



B. Flow structure and statistical properties of the velocity field 

PIV and LDV techniques have been used in the past in various experiments and setups 0, 0, HE ED] to 
characterize the structure of the flow field, to obtain velocity profiles and to infer and visualize the topology of the 
elastic von Karman flow. The visual characterization of the flow structure in a regime of elastic turbulence has been 
initially obtained from a few snapshots of the flow, by seeding it with light reflecting flakes (1% of Kalliroscope 
liquid) Q . Already such qualitative visualization revealed the presence of a big persistent toroidal vortex, which had 
the dimension of the whole setup. The average flow velocity along the radial direction was measured by LDV in a few 
points, and the results were quite well compatible with the presence of the big vortical structure 

With PIV we obtained a more complete and quantitative characterization of the global flow structure and of its 
dynamics. Some instantaneous vector fields of the horizontal components of the velocity, as well as the average of 
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Figure 6: (a) scaled PDFs of the injected power in constant SI forcing regime for different Wi. (b) scaled PDFs of the power 
fluctuations in constant T forcing regime for different Wi. The symbols (color online) are: down triangles (blue)-Wi = 5, up 
triangles {gie&a)-Wi = 17, circles (ved)-Wi = 25, diamonds (light blue)-Wz = 34, squares (black)-Wi = 46. 
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Figure 7: (a) Normalized spike periodicity, T s /A, in time series of the injected power vs. Wi: circles-constant SI mode, squares- 
constant T mode, (b) Normalized spike periodicity, r a /A, in time series of the azimuthal velocity vs. Wi. Data were collected 
at r = 2R c /3, z = d/2 in setup 2. 



2000 instantaneous fields taken for several values of Wi at a middle distance between the plates in the setup 1 are 
presented in Figs. El and The upper row of Fig.Elshows the instantaneous vector fields at three increasing values of 
Wi, while the lower row presents averaged vector field at the same Wi. One can easily identify the core of the toroidal 
vortex that appears in all images above the threshold of the instability, and a spiral vortex that additionally occurs 
in the elastic turbulence regime. These flow structures and their reorganization can be also observed in separate 
presentations of the averaged azimuthal and radial components, and rms of fluctuations of the azimuthal velocity 
component of the velocity field at the same Wi (see in Fig. I10I11JI . Below the elastic instability the profile of the 
average azimuthal velocity, Ve, exhibits a linear increase along the radius with a slope Q~ x which corresponds to a 
rigid body rotation (see the inset in Fig. 112b). while there is no motion in the radial direction, V r ps (Fig. 112b). 
Above the onset of the instability one can clearly see the creation of the core of the toroidal vortex at the cell center, 
and the restructuring of the radial motion [^3, ED] . 

Large toroidal vortexes driven by the hoop stress are actually quite well known to appear in swirling flows of elastic 
fluids pi 13 • The inhomogeneity of the shear rate profile in the primary laminar flow has long been recognized as 
their common origin. In our system this vortex first arises as a stationary structure at low shear rates, and leads to 
some increase in flow resistance even before the elastic instability 0|. Therefore, one concludes that the transition to 
elastic turbulence in the swirling flow between two plates is mediated by this vortex. The toroidal vortex is providing 
a smooth, large scale velocity field (see lower panels in Fig. EJ, which is randomly fluctuating in time (upper panels 
in Fig. and in which the fluid and the embedded stress tensor are chaotically advected. This type of advection 
can create variations of stress in a range of smaller scales, which may cause small scale fluid motion 0, E3- This is 
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Figure 8: (a) Normalized average spike width, c a „/A, vs Wi. (b) Normalized rms of the spike width, c rms /A, vs Wi. Data were 
collected in constant SI mode. 
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Figure 9: Instantaneous (upper row) and time averaged (lower row) velocity fields for different values of Wi: (a) Wi = 2.48, 
(b) Wi = 9.88, (c) Wi = 18.96. The data were collected in setup 1 at middle distance between plates. 



analogous to the generation of small scale variations of scalar concentration in chaotic mixing by large fluctuating 
vortexes. 

Above the instability threshold, the large scale toroidal vortex is also responsible for the ring-shaped topology of the 

— -1/2 

fluctuation field of the azimuthal velocity, V g rms = Vi (Fig. Illjl . The velocity fluctuations in the laminar regime 
are only due to instrumental noise (panel (a) in Fig. lTT|l . The average radial velocity changes its sign at about a half 
of the radius (see Fig. IT2"h). At higher values of Wi, the toroidal vortex forces the transition to elastic turbulence. 
This is accompanied by a second reorganization of the flow structure, with the appearance of a fluctuating spiral 
vortex (Figs. lflTlOjl and takes part in the scaling regime of P, P rms (see Fig. 01. Movie animations of the velocity fields 
evidence that this secondary fluctuating vortex is carried around and "rides" on top of the primary vortical flow. At 
the same time the circular symmetry of Vg ms is broken into a dipolar one (see Fig. lTT|l . The radial profiles of Vg and 
V r also change drastically. The structural transitions in the flow are also reflected in the average radial gradients of 
the azimuthal velocity and in the average vorticity field, as presented in Fig.^|for various values of Wi. 

Another way to characterize a turbulent flow is to display the radial profiles (averaged over the azimuthal coordinate) 
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Figure 10: Distribution of the radial velocity (left column) and the azimuthal velocity components (right column) for different 
values of Wi: (a) Wi = 2.48, (b) Wi = 9.88, (c) Wi = 18.96. The data were collected in setup 1 at middle distance between 
plates. 

□ □ □ 



C3 



xlO (m/s) 



m 



Figure 11: Fields of fluctuations of the azimuthal velocity at different Wi: (a) Wi = 8.32, (b) Wi = 9.88, (c) Wi = 11.1, (d) 
Wi = 12.72, (e) Wi = 13.83, (f) Wi = 17, (g) Wi = 19. The data were collected in the setup 1 at middle distance between 
plates. 

of the turbulent intensity, denned as It = V " Vfj at several values of Wi (see Fig. HOI- The velocity fluctuations 
in the laminar regime occur only due to instrumental errors. Above the elastic transition I t increases sharply but 
remains rather uniform on the level between 20 and 30% in a peripheral region for the radius ratios, r/Rc, between 
0.2 and 1. In the regime of elastic turbulence It further increa ses, and its dependence on Wi, presented in the inset 
in Fig.El exhibits the power-law scaling, I t - Wi 0A9±0 - 06 |3Sll3flj. 

The PIV measurement of time dependent velocity fields allowed us to calculate the average velocity gradients 
and vertical vorticity, and rms of their resjjective fluctuations, without involving the Taylor hypothesis that can be 
questionable for a smooth random flow 0, The typical radial distribution of the velocity gradients and of 

vorticity, averaged spatially over the azimuthal coordinate and temporally over 2000 images is rather uniform in the 
bulk for all Wi, but increases sharply near the wall (see Fig. Ilftt a)). while the rms of the velocity gradients and 
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Figure 12: (a) Profiles of the average azimuthal velocity, Vg, for several values of Wi. The inset shows a typical laminar profile 
of the azimuthal velocity component. (b)Profiles of the average radial velocity, V r , for several values of Wi. The labels are: 
(1)-Wi = 2.48, (2)-Wi = 4.41, (3)-Wi = 11.1, (4)-Wi = 15, (5)-Wi = 18. Data were collected in set-up 1 at middle distance 
between plates. 
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Figure 13: Average vertical component of the vorticity, u) z , at different Wi: (a) Wi = 8.32, (b) Wi = 9.88, (c) Wi = 11.1, (d) 
Wi — 12.72, (e) Wi = 13.83, (f) Wi = 17, (g) = 19. Data were collected in setup 1 at middle distance between plates. 
The gridded pattern visible in panels (a), (b) is an artifact of combined peak locking and numerical differentiation effects and 
should be disregarded. 



vorticity gradually increases with the radius (Fig, llfif a)). The dependence of the average vertical vorticity on Wi is 
displayed in Fig. HST h). At higher values of Wi one finds a gradual growth of the vorticity both in the bulk and near 
the wall. On the other hand, the plot in Fig.EJb) shows the Wi dependence of the rms of the vorticity, u) rms , scaled 
by A, at several locations along a radius in the bulk. The scaled rms of the vorticity saturates in the elastic turbulence 
regime at a Wi between 12 and 35, that is consistent with the recent theoretical prediction. Similar behavior is 
shown by all the components of the velocity gradient, and is consistent with the theoretical predictions, though the 
saturation value uo r z ms ■ A ~ 2 is higher than expected @, 0] ■ The latter means probably that the nonlinearity of the 
polymer elasticity also contributes to the saturation of the elastic stress jUI^] (see further discussion of this issue). 

Time series of azimuthal velocity obtained from LDV measurements for several values of Wi are displayed in 
Fig- EI a). As Wi increases, the velocity signal becomes increasingly intermittent. PDFs of the normalized azimuthal 
velocity, (Vg — V g av )/V g rms , are shown in Fig. IT7T hV The right side skewness of the PDFs is due to the rare events 
(spikes) which intermittently occur in the velocity signal (Fig. I17f a)). which are probably related to the dynamics of 
the randomly fluctuating spiral vortex previously discussed. 

In order to understand better the structure of the flow in the regime of elastic turbulence, we have studied the 
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Figure 14: Turbulence intensity versus the reduced radial coordinate at different Wi: diamonds- Wi = 31.56, squares- Wi = 
27.73, circles- Wi = 23.67, up triangles-Wi = 19.29, down triangles-Wi = 15.48, stars-Wi = 7.54, half filled triangles- = 2.82. 
The inset shows the dependence of the turbulence intensity on Wi at r/R c — 0.85 obtained from LDV measurements. The full 
line is a power law fit, It oc Wi°' 49±0 ' 0(i . The arrow indicates the onset of the elastic instability. The data were collected in 
setup 2. 
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Figure 15: (a) Profiles of the vorticity (lj z ) averaged spatially in the azimuthal direction and temporally over 2000 images, for 
different Wi. The labels are: (l)-Wi = 2.82, (2)-Wi = 5.4, (3)-Wi = 10, (4)-Wi = 14.48, (5)-Wi = 17.43, (6)-Wi = 19.29, 
(7)-Wi = 26.13, (8)-Wi = 27.73, (9)-Wi = 32.3, {10)-Wi = 34.5. (b) Dependence of the average vorticity on Wi at 
different radial positions r/R c ,: left triangles-0.1, down triangles-0.2, squares-0.33, up triangles-0.4, stars-0.5, circles-0.66, right 
triangles-0.7, half filled squares-0.8, half filled circles-0.9. The arrow marks the onset of the elastic instability. The data were 
collected in setup 2. 



dependence of the statistics of the azimuthal velocity on the radial coordinate, for a fixed value of Wi = 36.1 and a 
given vertical coordinate, z = d/2. As illustrated in Fig. ITsT al. in the central (r = 0) and peripheral (r = 15i? c /16) 
flow regions, the velocity signal shows no significant intermittency, whereas around r = R c /2 it is strongly intermittent. 
The radial dependence of the level of intermittency reflects the presence of the randomly fluctuating spiral vortex 
(where, probably, r = R c /2 corresponds to the arm of the spiral). This behavior is also reflected in the PDFs of the 
normalized azimuthal velocity presented in Fig. HsT b): near the center and the vertical wall of the cell, the distributions 
are symmetric and single peaked whereas around r = R c /2 they become strongly skewed and doubly peaked. 

Fig. 1191 shows PDFs of the normalized accelerations as calculated from the temporal LDV measurements, (Y(t) — 
Y av )/Y rms , where Y(t) = dV d f' V^ 1 . We note that all the data for several values of Wi collapse on a single curve 
remarkably well, and that moreover the PDFs show clear exponential tails. 

The PIV data allow us also to calculate various correlation and structure functions of velocity, velocity gradients, 
and vorticity. Their scalings can give information on the degree of deviation from a Gaussian random field, a standard 
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Figure 16: (a) Profiles of the rms of the fluctuations of the vorticity, u™ 8 , at different Wi. The labels are: (l)-Wi = 2.82, 
(2)-Wi = 5.4, (S)-Wi = 10, (4)-Wi = 14.8, (5)-Wi = 17.43, (6)-Wi = 19.29, (7)-Wi = 26.3, (8)-Wi = 26.3, (9)-Wi = 27.73, 
(lO)-Wi = 32.3, (b) Dependence of the scaled rms of the vorticity, uj rm3 ■ A, on Wi at different radial positions r/R c : full 
squares, 0.2; open squares, 0.33; full circles, 0.4; open circles, 0.5; diamonds, 0.66. The data were collected in setup 2. 
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Figure 17: (a) Time series of azimuthal velocities, Vg, for several values of Wi. (b) PDFs of normalized azimuthal velocities 
for various Wi: squares-19.6, circles-21.7, up triangles-24.5, down triangles-27, diamonds-29.6, left triangles-33.2. Data were 
collected in setup 2, at z = d/2 and r = 2R c /3. 



way in hydrodynamics to quantify intermittency and to compare it with the corresponding scaling in known cases. 

The spatial power spectra of the velocity were obtained from PIV data by averaging 2000 instantaneous spatial 
spectra. Although the finite spatial resolution of PIV, which limits the accessible range of wave numbers, leads to 
an artificial cut-off at k ~ 3000m -1 , a power law decay k~ s with 8 ~ 3.5 is clearly observed (Fig. 1201 [3^. Igg lh This 
decay is consistent with the scaling obtained earlier for the velocity power spectra in the frequency domain 
The large value of the exponent 8 implies that the power of fluctuations decays very quickly as the size of eddies 
decreases. The main contribution to the fluctuations of both the velocity and the velocity gradients is therefore due 
to the largest eddies, and the power of the latter should scale as fc - '* 5-2 '. 

The agreement between the velocity spectra measured at a single point in frequency domain 0, 0| and the spectra 
measured directly in fc-domain pSfll. l33| deserves a brief discussion. Following Lumley ^3], the relation between the 
spatial spectrum P{k) and the frequency spectrum P(f) can be written as 



I? d 2 (k 2 P(k)) 

p( k ) = v.p(f) - -± ■ K dk ; ]) +o(it) , 



(8) 
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Figure 18: (a) Time series of azimuthal velocities, Vg, for Wi = 36.1 at several radial positions, r. (b) PDFs of the normal- 
ized azimuthal velocity, for Wi = 36.1, at several radial positions r: squares-0,circles-i? c /8,up triangles- R c /2, down triangles- 
37? c /4, diamonds- 7R c /8, left triangles-15i? c /16. The data were collected in setup 2, at z = d/2. 
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Figure 19: PDFs of normalized accelerations for various Wi: full circles-14; open circles-16.2; full squares-19.8; open squares- 
22.3; up-triangles-26.3; down-triangles-30.7; diamonds-38.9. The data were collected in setup 2 at z = d/2 and r = 2R c /3. 
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where I t = Vg" ms /Ve, Vg" ms = Vg , and Vg are the rms of fluctuations of the azimuthal velocity and the average 
azimuthal velocity, respectively. If P(k) oc k 52 and P(f) oc f Sl , the equation above leads to: 

1 + f • 5 2 ■ (S 2 + 2) 



log 
5i oc 



log(fc) 



(9) 



If one plugs into the last equation 62 ~ —3.5, the difference between the exponents is (for k w 1000m -1 ) as small as 
$2 — Si r* 0.1. Thus, the experimental resolution does not allow us to observe the difference in the scaling exponents 
for the spatial and temporal spectra 

Next, let us focus on the typical spatial and temporal correlation scales in elastic turbulence. A typical scale, at 
which elastic stress is correlated can be estimated as L = 2ir J P(k)dk/ J k P(k)dk. In our setup and in the elastic 
turbulence regime, one gets L ~ 5.9 cm, i. e. on the order of the cell size. The Eulerian correlation time is defined 
as t c = / tC(t)dt/ J C(t)dt, where C(t) = V(T)V(T + t)/ (V rms f is the Eulerian correlation function. Figure EU 
presents the normalized t c /A as a function of Wi at some radial positions in the cell. The correlation time drops 
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Figure 20: Power spectra of the fluctuations of the azimuthal velocity component at different Wi: circles-Wi = 4.41, up 
triangles-M^i = 5.72, down triangles-Wi = 8.32, left triangles Wi = 11.1, right triangles-Wi = 12.7, hexagons- = 13.8, 
diamonds-Wi = 16, half filled squares- Wi = 18, empty cixcles-Wi = 19. The solid line is the power law AT 3 ' 5 . The data were 
collected in setup 1. 
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Figure 21: Scaled Eulerian correlation times of the azimuthal velocity component as function of Wi: squares-setup 1, circles- 
setup 2. 



significantly in the transition region and then saturates in the elastic turbulence regime, as the rms of the vorticity 
(and velocity gradients) fluctuations do. The related saturation of both the Lyapunov exponents of the Lagrangian 
trajectories, already reported (^J, and the reduced velocity correlation time, r c /A, shown in Fig. |2] (in the same 
range of Wi) just reinforces the observation. We point out that the inverse saturated value of correlation time is of 
the order of the saturated value of uj" ns and A" 1 , as predicted by theory 

Using the PIV data one can also calculate the structure functions of the velocity gradients, defined as 



S P (r) 



dVe(r a + r) dV e (r ) 



Or 



dr 



(10) 



and shown in Fig. |23 As for inertial turbulence, we can look for a scaling in the spatial range corresponding to the 
power law decay of the velocity spectra fFig. l2()Jl in the form: S p (r) = r*> p . The dependence of the normalized scaling 
exponent, (p/d, on the order of the structure functions, p, is presented in the inset of Fig.[22]for all flow parameters 
mentioned above. For comparison, the structure functions of the vorticity, of the injected power at constant fi, Pn, 
and constant torque, Pt, were also calculated, and the results of the analysis for the normalized scaling exponent 
were also plotted in the inset of Fig. |22j All these dependencies are rather close to that obtained for a passive scalar 
in elastic turbulence, presented on the same plot and based on our published data all of them deviate from 
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Figure 22: Structure functions of the velocity gradients of order p up to 6 as a function of the scaled radius. Inset: normalized 
scaling exponent vs. p for: injected power P«- full circles; injected power Pt- open circles; passive scalar- full squares; vorticity- 
open squares; velocity gradient- up-triangles. 



a linear dependence that is characteristic of a Gaussian distribution [45j . 

It is worth to make here some remarks about the applicability of the Taylor hypothesis to flows which are spat ially 
smooth random in time, such as is elastic turbulence, which we have studied experimentally in more detail in 
On one hand, the comparison of the power spectra of the velocity in frequency and wave number domain suggests 
that the Taylor hypothesis is applicable at least in regions close to the wall. On the other hand, from Fig.^Jone can 
learn that in the central flow region, at r/R c < 0.2, the rms of the velocity fluctuations becomes much larger that 
the mean velocity, fact which is normally thought to violate the assumptions of the hypothesis. By detailed analysis 
of cross-correlations as well as of the structure functions of the velocity fluctuations, it was indeed shown in that 
strong velocity fluctuations near the cell center prevent, and flow smoothness and lack of scale separations close to 
the boundaries limit the quantitative applicability of the Taylor hypothesis. However, this deficiency can be corrected 
by a proper choice of the advection velocity, so that a corrected Taylor hypothesis is applicable to most parts of this 
chaotic flow. 



C. Velocity and velocity gradient profiles and the boundary layer problem 

To learn more about the flow field generated by elastic turbulence, we investigated in detail the structure of 
horizontal and vertical velocity boundary layers for different values of Wi and fluid viscosity. The vertical profile of 
the azimuthal velocity, Vg , in the laminar flow of the pure solvent is nearly a straight line 0, H3| ■ The transition 
to elastic turbulence flow significantly changes the distribution of V e av , producing a high shear layer near the upper 
plate and a low shear region near the middle of the gap (at z/d = 0.5). Such a distribution of V g av is reminiscent of 
average velocity profiles in usual high Re turbulence. In Fig. |22l we present the average azimuthal velocity profiles 
measured by LDV in a vertical plane at r — 3i? c /4 of the setup 2 with d = 2cm height in 65% of saccharose solution 
and normalized by the velocity value at the disk, Vg{z/d) = V B av {z/d)/V g up . In Fig. I24I the same data for different Q 
are collapsed on a single curve, after subtraction the z-dependence of the average velocity in the bulk for each curve, 

that is V e *(z/d) = Vg(z/d) — Vg{\/2) — | d ^f/^ ^ , where d ^f/ff is obtained from a linear fit of the data points. 

Similar velocity profiles were also obtained with polymer solutions with saccharose content reduced to 60% and 50%. 
We defined the boundary layer width w as the intersection abscissa of the linear fits for the bulk and near-the-wall 
parts of the velocity profile. It is clear from the plot that w is independent of ft for each solution and is seen to depend 
only on the solution viscosity, whereas we obtained the fit w oc jy - 26 * 05 (upper inset of Fig. l24"jl . Fluctuations of the 
azimuthal velocity are small near the upper plate, reach a maximum at z/d ~ 0.1, and start to decrease at larger z. 
Again, such distribution of VJ ms along z is reminiscent of the velocity fluctuations in the boundary layer of turbulent 
flows of a Newtonian fluid ja, [2fj|. An example of the profile of the rms of the azimuthal velocity V e rms , obtained at 
various Vl in the setup 2 with height d = 1cm, is shown in Fig. [2^1 Well-defined peaks appear at the same position, 
which does not depend on ft, and match the boundary layer thickness w found with the previous analysis of the 



24 



average velocity profiles, and displayed in Fig, 1241 
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Figure 23: Vertical profile of the average azimuthal velocity, normalized by its maximum value at the disk, measured by LDV 
in setup 2, with d = 2cm, in a vertical plane at r = 3i? c /4, in a 65% saccharose polymer solution, for several values of fi: 
full squares (black)-5s _1 , circles (red)-4s _1 , up-triangles (green)-3s _1 , down-triangles (blue)-2.5s _1 , diamonds (cyan)-2s _1 , 
left-triangles (magenta)- 1 . 5s ~ 1 . 
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Figure 24: Compensated vertical profile of the average azimuthal velocity V e *(see text), for several values of fi: full squares-5s _1 ; 
circles-4s _1 ; up-triangles-3s _1 ; down-triangles-2.5s _1 ; diamonds-2s _1 ; open squares-1.5s _1 . Inset: boundary layer thickness w 
vs. solution viscosity; the solid line is the fit w oc }j - 26 ± 05 . Data measured in the same setup as data of Fig. I23I 

We also measured by PIV radial profiles of both the azimuthal and the radial velocity in a horizontal plane at 
z = d/2 in the vicinity of the vertical wall of the cell, shown in Fig. I26f a.b). As shown in Fig. I2fifa. the azimuthal 
velocity component Vg can be fitted by a linear function in the vicinity of the wall. This agrees with the requirements of 
incompressibility, boundary condition for the velocity and smoothness of the velocity field, as pointed out in Ref. 

In the regime of elastic turbulence the slopes of the linear fit functions increase monotonically with Wi (see the 
inset in Fig. I26b). Using the fits one can determine the width of the vertical (close to the exterior cylindrical wall) 
boundary layer, wg, of the azimuthal velocity. 

As shown in Fig. [23 the width of the vertical boundary layer, we, decreases abruptly above the onset of elastic 
instability and saturates in the regime of elastic turbulence, for large values of Wi. Unfortunately the limited accuracy 
of measurements of the radial velocity component near the wall (note that radial velocities are about two orders of 
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Figure 25: Profile of the rms of the azimuthal velocity measured by LDV in a vertical plane at r = 3i? c /4, in setup 2 with 
d = 1cm height, in 65% of saccharose solution, for several values of f2: full squares-6r-ad/s, full circles-5.5rod/s, full up triangles- 
5rad/s, full down triangles-5.5rad/s, empty squares-4rad/s, empty circles-3.7rad/s, empty up triangles-3.5rad/s, empty down 
triangles- 3rad/s, full hexagons-2.5rad/s, empty hexagons-2rad/s. 
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Figure 26: Near-wall profiles of the velocity in a horizontal plane at z = d/2, for various values of Wi in setup 2, with 65% 
saccharose polymer solution, (a) Radial profiles of the tangential velocity component V$. The full lines are linear fit functions. 
The inset shows the dependence of the slope of the linear fit functions, s, on Wi. (b) Radial profiles of the radial velocity 
component V r . The symbols are: full squares-Wi = 2.4, full up triangles-Wi = 4, full diamonds-Wi = 5.9, full right triangles- 
Wi = 7, full stars-Wi = 10.6, empty circles-Wi = 16.6, empty squares-Wi = 20.3, empty up triangles-Wi = 23.8, empty 
diamonds- Wi = 23.8, empty right triangles- Wi — 30.2. 



magnitude smaller than the azimuthal ones) did not allow us to assess beyond doubt the functional behavior of the 
radial profiles of radial velocity component near wall, presented in Fig. 126b . 

In order to improve the spatial resolution of the velocity field measurements near the wall, and particularly to resolve 
the horizontal boundary layer by looking at the rms of the radial gradient of the azimuthal velocity, (dVe / dr) rms , 
the field of view of the PIV measurements was reduced down to 10mm in the radial direction. The radial profiles 
of {dVo/dr)™rm near the wall for various values of Wi in the elastic turbulence regime are shown in Fig. 1281 As in 
the case of the average velocity in the vertical plane (Fig. 12411 . the bulk part of the radial dependence fitted by a line 
was subtracted for each curve. Well-defined peaks, positioned independently of Wi, can easily identified inside the 
velocity boundary layer. We would like to point out that the value of (dVe/dr) rms at the maximum exceeds up to 
two orders of magnitude its bulk value (see Fig. Ilfif b)). 

The same data were also used to determine the average value of {dVg / dr) rms in the boundary layer. The latter was 
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Figure 27: Dependence of the width of the vertical boundary layer, wg (defined via the width of the linear part of the near wall 
PIV profiles of the tangential velocity component) on Wi. Data were collected in setup 2. 
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Figure 28: Normalized radial profiles of rms of the azimuthal velocity gradients, {dVg/dr)™'m in the vicinity of the vertical 
wall, for various values of Wi (color online): full squares- Wi = 10.6, empty squares- Wi = 12.4, full up triangles- Wi = 16.5, 
empty up triangles-Wi = 18.1, full circles-Wi = 20.3, full left triangles-Wi = 22, empty circles-Wi = 23.8, empty left triangles- 
Wi = 25.4, full right triangles-Wi = 27, empty right triangles-Wi = 28.6, stars-Wi = 30.2. Data were collected using PIV in 
setup 2 with a 65% saccharose polymer solution, in a horizontal plane at z = d/2. 



obtained from widths of distributions of the slopes of the radial dependence of the azimuthal velocity in the boundary 
layer region, for each value of Wi. These average values of (dVg/dr) rms scaled by A, in the vertical boundary layer of 
the setup 2 with cell height d = 1cm are shown in Fig. [221 as function of Wi. A clear increase of {^-) rms • A 3> 1 in 
the boundary layer, in contrast to its saturation value in the bulk, {^-) rms ■ A « 2, was found in the same range of 
Wi (see further discussion). 



D. Discussion: role of the elastic stresses in elastic turbulence in a finite size cell. 



Let us summarize our different experimental observations, discuss some of their implications, and delineate a 
comprehensive picture of the phenomenon called elastic turbulence. 
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Figure 29: Scaled rms of the velocity gradients in the boundary layer vs Wi. The arrow indicates the onset of the primary 
elastic instability. 



(i) The rms of the vorticity (and that of velocity gradients) saturates in the bulk of the elastic turbulent flow, 
leading to the saturation of the elastic stress. The saturation of u) rms at high values of Wi in the elastic turbulence 
regime, observed in the data in Fig. gives us the possibility to test which of the two theoretical mechanisms is 
responsible for it: either the feedback reaction of polymer molecules on the flow [24[ or the nonlinearity of polymer 
elasticity The local Weissenberg number can be estimated from Fig. EH as Wi 1 = u rma ■ A ~ 2. which exceeds 
the unit value predicted by the model with feedback of linear molecule elasticity, discussed by 0,113] • We conclude 
that the nonlinearity of the polymer elasticity may also contribute to the saturation of the elastic stresses pS]. The 
saturation of the elastic stress in the bulk also naturally explains the power-law behavior of the average injected 
power, -p^— ~ Wi 0A9 , seen in the plot in Fig.GD Indeed, the injected power is proportional to the torque, P = T fl. 
The latter is just the shear stress averaged over the upper plate, and its growth above the instability threshold occurs 
solely due to the elastic stress, t p , @. In the elastic turbulence regime, however, the theory predicts t p ~ rj/X. On 
the other hand, it was found that, due to shear thinning of the solution used, A ~ fl~ - 3 . Thus, the theoretically 
expected growth of the mean injected power with respect to its laminar value, Pi am , should be solely due to the elastic 
stress and have the following power-law scaling: ~ t p ~ Wi 0A3 which is rather close to the scaling we observed 
experimentally. 

Next, we made three main observations (ii)-(iv) below, which extend further the formal similarity with hydro- 
dynamic turbulence, despite the distinctively different underlying mechanism for intermittency. Our findings in the 
statistics and structure of the velocity field have to be related to the finite size of the experimental vessel, whereas 
the current theory describes elastic turbulence in an unbounded domain 0, 0] . One can expect to find new effects in 
statistical and scaling behavior, similar to those recently studied both theoretically and experimentally for a passive 
scalar advected by a smooth random flows in bounded geometries (Ts j. such as that in a channel considered by (T3 |. 

As was already pointed out in the Introduction, there is a deep analogy in the statistical behavior of elastic 
turbulence, of the fast magneto-dynamo, and for the passive scalar dispersion in the Batchelor regime. As was shown 
in Section III Dl the elastically driven random flow can be considered as turbulence of the elastic stress field; all its 
statistical properties can be explained by the interplay of (a) stretching and folding of fluid elements carrying elastic 
stresses on Lagrangian trajectories in the Batchelor regime, (b) dynamo effect leading to the saturation of stresses 
in the bulk flow, and (c) homogeneous linear damping of the stress. Thus, since the dynamics of the elastic stress 
tensor can be reduced to the dynamics of the vector field B, this scenario is very close to that found in the fast 
magneto-dynamo except that (c) is replaced by the diffusive dissipation of the magnetic field. On the other hand, 
the passive scalar problem differs from elastic turbulence (besides the fact that in this case a vector is replaced by a 
scalar) in the item (b), since in the former the feedback reaction on the flow is absent. However, in spite of the fact 
that there is no exact similarity between these three problems, the main ingredient (a) gives a basis for the similarity 
in statistical and scaling behavior of these systems. Taking into account the similarity, we can extend it to statistical 
properties of elastic turbulence in a finite size vessel. 

/ AT/ \ VTYLS 

(ii) The rms of the velocity gradients y^jr) or oj rms , and thus the elastic stress, grows linearly with Wi in 
the boundary layer (see Fig. l29j) . near the driving disk. The rms of the velocity gradients saturates in the bulk (see 
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Fig. ll6|l , Moreover, the values in the boundary layer are one to two orders of magnitude larger than in the bulk. All 
this implies that a boundary layer is formed, with an excess of the elastic stress is concentrated in it. 

(iii) The PDFs of the injected power SPq/P™ 3 , at either constant angular speed or torque, show skewness and 
exponential tails, which both indicate intermittent statistical behavior. Also the PDFs of the normalized accelera- 
tions, which can be related to the statistics of velocity gradients via the Taylor hypothesis, exhibit well pronounced 
exponential tails. We can match this facts and the emergence of the boundary layer with the following scenario: the 
forcing from the moving boundary (the upper disk) produces an accumulation of the elastic stress in the boundary 
layer due to a constant flux of momentum from the upper boundary. Elastic stress is intermittently injected into 
the bulk by random release events. Whenever that happens, the power on the disk decreases. Constant build-up 
and sudden release would thus explain the skewness of the PDF toward small values of the injected power, while 
intermittency would be the origin of the observed exponential tails. The exponential tails in the PDF of the velocity 
gradients and accelerations (Fig. 1 1 71 1 911 can be also explained in the same way, except that for them the symmetry of 
the PDF is preserved due to injections from both the upper and lower boundary layers. 

However, the most surprising and striking feature of elastic turbulence observed is (iv), the presence of a velocity 
boundary layer and the emergence of a new length scale. Since elastic turbulence is a turbulence of an elastic stress 
field, one can expect primarily the appearance of a boundary layer in the elastic stress field. We observe in addition 
the presence of a velocity boundary layer, which just reflects some of the features of the boundary layer of the elastic 
stresses (rms of the velocity gradients) . The measured thickness of the velocity boundary layer is much smaller than 
the vessel size. To argument further on it, we note that the picture presented in point (iii) is in close analogy with 
the turbulent advection of a passive scalar in the Batchelor regime in a finite channel flow, where the excess of tracer 
from the boundary layer is intermittently injected into the bulk 0,0]. There the boundary layer for a passive scalar 
(mixing boundary layer) results from zero boundary conditions for the flow velocity components along and towards 
the wall. This leads to less effective advection of the passive scalar near the wall, causing an accumulation in the 
near-wall region, and thus a boundary layer in which the statistics of the scalar is different than that in the bulk flow. 
In the case of the passive scalar, the mixing boundary layer and the new length scale appear due to the existence 
of a small parameter, the inverse Peclet number, Pe _1 = (D/V rms L) <C 1, where D is the diffusion constant, and 
L is the vessel size. Then, as predicted theoretically an d confirmed experimentally the mixing boundary 
layer width is found to scale as w oc Pe~ 1//4 oc D 1 / 4 . Intermittent ejection of the tracer into the bulk, analogous to 
the intermittent injection of elastic stresses in elastic turbulence, reduces the effectiveness of mixing, increases the 
mixing length, and causes stronger (power-law) dependence of the mixing length on Pe than in an unbounded vessel 
(in which dependence should be logarithmic). 

Carrying on the analogy with the passive scalar problem, one can suggest that also here a small parameter of the 
problem defines a new length scale. Since diffusion does not play any significant role in elastic turbulence, we suggest 
that 



V dr /bulk/ V dr J bl 

can be considered as the relevant small parameter. This non-uniform distribution of the elastic stresses should result 
in, and be directly related to, a non-uniform distribution of polymer stretching: the polymer molecules should be 
stretched considerably more in the boundary layer than in the bulk. We have verified this prediction in a separate 
experiment on a single molecule in a random flow Thus, the width of the boundary layer for the elastic stress, 
based on these considerations and on the experimental observations, scales like w/L oc N^r] 1 / 4 , with \x an exponent 
to be determined. While it is appealing to note that the exponent of r\ (Fig. 12} is the same as that of 1/Pe in 
the previous case, we are not really able to justify it theoretically. On one hand, in contrast with the passive scalar 
problem, where dissipation due to diffusion introduces the natural length scale, elastic stresses in elastic turbulence 
are dissipated via linear damping with the polymer relaxation time, A. Since A cx rj, r\ has good reason to appear in 
the scaling relation for w, but we cannot guess a value of its exponent. On the other hand, the suggested scaling has a 
serious flaw: while we have seen that (^r) b lk saturates in the regime of well developed elastic turbulence (Fig.H?)h). 

the suggested scaling law bears a dependence on fl via (^f ) b \ (see Fig. I29fl . whereas we found experimentally that 
w is independent of it. These two issues need to be addressed by a forthcoming theory. 

Finally, (v), the scaling of the structure functions of the vorticity, velocity gradients, and injected power is found 
to be the same of that of a passive scalar advected by an elastic turbulent velocity field (see subsection IIV BJl , fact 
which further reinforces the analogy in the statistics and scaling of the passive scalar and of the elastic stresses due 
to advection by a smooth random flow. 
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V. CONCLUSIONS 

The main message of this study is a surprising similarity in scaling, statistics, and structure of both the elastic 
stresses and the passive scalar in elastic turbulence in a finite size vessel, in spite of the important difference in the 
dynamo effect. 

First of all, we have confirmed that the flow is smooth and essentially dominated by large scale structures: the 
velocity power spectra decays asP~ A: -3 3 , which suggests that the power of fluctuations of velocity gradients scales 
as fc~ 13 (see Sec. HVB)) . As explained in Sec. Ill Dl the main contribution to the fluctuations of the velocity gradients 
and of the velocity differences at all scales comes from the largest eddies, having dimensions of the whole system (the 
gap between plates). This conclusion has an immediate implication for the mixing in the flow: it should result in the 
same type of patterns and in the same statistics observed in a randomly varying in time flow with a linear velocity 
field (locally uniform rate of strain), V(r, t) = Vn(£) + • (r — r(0)), which approximates well the velocity field in 
elastic turbulence 0] (r is the position vector). 

Looking at more detail into the local properties of the velocity field, we have evidenced a saturation of the rms 
of the vorticity (or velocity gradients), which should lead to the saturation of the elastic stresses in the bulk. This 
agrees with theoretical predictions, save the fact that the saturation level found is about twice than what theoretically 
predicted The discrepancy can be attributed to nonlinear polymer elasticity |23 |. The observed saturation is 

however remarkably consistent with the scaling of the averaged injected power with Wi. 

The presence of a boundary layer for stress and velocity was also evidenced and characterized: the rms of the 
velocity gradients (and therefore the elastic stresses) grow linearly with Wi within it, and exceed the corresponding 
value in the bulk by one-two orders of magnitude, suggesting that the elastic stresses are accumulated near the wall 
and are intermittently injected into the bulk. The view is confirmed by the properties of the PDFs of the injected 
power and by the exponential tails of PDFs of the normalized accelerations. This leads to a model for the dynamics 
of elastic turbulence, in which elastic stress is introduced into the fluid by the driving boundary, accumulates in the 
boundary layer and is intermittently injected into the bulk of the flow. The situation is entirely similar to the trapping 
of a passive scalar in the mixing boundary layer in a finite size vessel, and to its intermittent injection into the bulk. 
The non-uniform distribution of the elastic stresses, which causes non-uniform polymer stretching is currently verified 
experimentally (-L^ . 

Universality and scaling of the boundary layer thickness are also investigated. The emergence of a relevant scale for 
it, based on properties measured from the velocity profile, is the most unexpected and striking result, since according 
to the theory of elastic turbulence in unbounded domains, there would be no other length scale besides the dissipation 
length or the vessel size. The measured boundary layer is much thinner than the vessel size, and the scaling of its 
width w with 77 was found experimentally. A small parameter, which is related to the polymer extension ratio in the 
bulk and in the boundary layer, was identified in this problem. 

Finally, the similarity between elastic turbulence and Batchelor regime advection of passive scalar is reinforced by 
the scalings of structure functions of the vorticity, velocity gradients, and injected power, which were found to be the 
same as that for a passive scalar advected in a velocity field produced by elastic turbulence. 

The experimental results presented call for further development in the theory of elastic turbulence in a container of 
finite size. Particular attention should be devoted to what causes the scaling of the boundary layer width, the scaling 
of structure functions and rms of the injected power, and the saturation level of the elastic stresses in the bulk flow. 
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